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Introduction: Combustion instability, where 

unsteady heat release couples with acoustic modes, has 
long been an area of concern in liquid rocket engines. 
Accurate modeling of the acoustic normal modes of the 
combustion chamber is important to understanding and 
preventing combustion instability. The injector resistance 
can have a significant influence on the chamber normal 
mode shape, and hence on the system stability. 
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Figure 1. System geometry 


Computational Methods: This study evaluates 

the effect of injector resistance on the mode shapes and 
complex eigenfrequencies of an injector/combustion 
chamber system by defining a high Mach-flow form of 
the convective wave equation (see Eq. 1) in COMSOL 
Multiphysic’s Coefficient Form PDE Mathematics 
Module. The form of the wave equation shown in Eq. 1 
is based off of a similar form derived by Campos in Ref. 
[1]. The additional terms in Eq. 1 are included by adding 
source terms to COMSOL’s base governing equation. 
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Figure 2. Pressure drop normalized by manifold pressure 
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Figure 3. Axial velocity, normalized by maximum 

Results: A larger pressure drop/velocity across 
the orifice leads to higher acoustical resistance. 
Although there is only a 5% difference in pressure 
drop between the two cases, there is a significant 
impact on the acoustic mode shape of the 
system. For the low resistance case, the 
manifold/injector boundary behaves as an 
acoustically open boundary, while for the high 
resistance case it behaves as acoustically closed. 
Additionally, increasing the resistance across the 
injector leads to a change in both the real and 
imaginary components of the complex 
eigenfrequencies, as shown in Figure 4. 
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c 0 = Speed of sound 

p 0 = Density 

v 0 = Velocity vector 

A = a + ico, Complex eigenvalue 

m 

ip = ip r + iip 1 , Complex velocity potential 

Background steady-state flow conditions are 
determined through a NASA Marshall Space Flight 
Center in-house computational fluid dynamics model, 
and interpolated onto the COMSOL mesh. Two cases 
are investigated, one with an injector orifice diameter 
of 6.502 mm (0.256 in), and the other with a 
diameter of 5.613 mm 0.221 in). As shown in Figure 
2, the larger orifice leads to a 28% drop in pressure 
with respect to manifold pressure, while the smaller 
orifice provides a 33% drop. 


A. Large orifice, injector B. Small orifice, injector 

resonance, 25 + 9219i resonance, -395 + 6904i 



C. Large orifice, chamber D. Small orifice, chamber 
resonance, 106 + 14411i resonance, -70 + 12547i 


Figure 4. Sample mode shapes for the small and large orifices 

Conclusions: The knowledge gained through 

this model can be used during future design cycles to 
favorably shape the combustion chamber mode 
shape, and to determine the complex 
eigenfrequencies in an effort to predict which modes 
are susceptible to instability. 
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